Entanglement in a second order quantum phase transition 
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We consider a system of mutually interacting spin 1/2 embedded in a transverse magnetic field 
which undergo a second order quantum phase transition. We analyze the entanglement properties 
and the spin squeezing of the ground state and show that, contrarily to the one-dimensional case, a 
cusp- like singularity appears at the critical point A c , in the thermodynamical limit. We also show 
that there exists a value Ao > A c above which the ground state is not spin squeezed despite a 
nonvanishing concurrence. 
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Entanglement is a truly specific property of the quan- 
tum world, and one of its deepest signature, as it was al- 
ready recognized in the early days of quantum mechanics. 
It is at the heart of the celebrated EPR paradox 0, Q , 
and plays a key role in the measurement problem as well 
as in the quantum to classical transition Entangle- 
ment is also central in quantum computation 0,0 where 
the most interesting operations cannot be completely ful- 
filled through the manipulation of separable states. 

Recently, entanglement properties of systems undergo- 
ing quantum phase transitions Q have attracted much 

attention[lH|l[ll[l3. 

Interestingly, the concurrence 
of the ground state which is related to the entanglement 
of formation ^3] , has been shown to be strongly affected 
at the critical point 

Hi 

More precisely, in the one- 
dimensional (ID) Ising model in a transverse field, Os- 
terloh et al. have shown that the derivative of the concur- 
rence with respect to the coupling constant diverges at 
the transition point jg| although the concurrence itself is 
not maximum. These pionneering results raise the ques- 
tion of the universality of these behaviors. Apart from 
ID quantum spin models, there has been, up to now, no 
other analysis of the ground state entanglement in sys- 
tems displaying quantum phase transitions except in the 
Kagome lattice |9| . Actually, the lack of exact solutions 
especially in higher dimensions implies a numerical treat- 
ment which often restrict the study to a small number of 
degrees of freedom. Such approaches do not allow, in 
general, an accurate description of the thermodynamical 
properties. 

In this Letter, we study the entanglement properties 
of a quantum system made up of N spins 1/2 on a sim- 
plex (each spin interacts with all others) embedded in a 
magnetic field. The permutation symmetry of this sys- 
tem allows us to restrict the ground state determination 
to a iV-dimensional subspace and, consequently, to deal 
with a large number of spins (about one thousand). We 
analyze the concurrence and the spin squeezing of the 
ground state which are, in this case, closely related [l3j . 
Contrarily to what happens in the ID Ising model, the 
concurrence of the ground state is maximum and displays 
a cusp-like singularity at the critical point. Moreover, at 



the transition point, the ground state is maximally spin 
squeezed and its squeezing parameter 0] vanishes in the 
thermodynamical limit. Finally, we show that there ex- 
ists a special line in the parameter space where the con- 
currence vanishes and above which the ground state is 
not spin squeezed although the concurrence is nonzero. 

Let us consider the following Hamiltonian first intro- 
duced by Lipkin et al. |lf>l |: 

h = -4e(«+t«)-e<' o) 



i<j 



= _^(S2 + 7 S«)-2S, + ^(1+ 7 ) I (2) 

where the cr Q 's are the Pauli matrices and S a = J^i a a/2- 
We focus here on the ferromagnetic case (A > 0) and we 
mainly consider the case < 7 < 1. The prefactor 1/N 
is necessary to get a finite free energy per spin in the 
thermodynamical limit. 

The Hamiltonian H preserves the total spin and does 
not couple states having a different parity of the number 
of spin pointing in the magnetic field direction, namely: 



[H,S 2 ] = 0, 



= 0. 



(3) 
(4) 



for all 7. In the isotropic case 7=1, one further has 
[H, S z ] = so that H is diagonal in the standard eigen- 
basis {1-5, M)} of S 2 and S z . 

For any 7, this system displays a second order quantum 
phase transition at A c = 1 which is characterized by the 
mean- field exponents 0, ^3 ■ Nevertheless, a mean-field 
approach cannot provide nontrivial entanglement prop- 
erties since it essentially turns the Hamiltonian into a 
sum of single-body Hamiltonians. It is thus necessary to 
use numerical diagonalizations of H for finite N. The 
dimension of the Hilbert space is 2 N but the study of the 
ground state reduces to a problem linear with N since it 
lies in the fully symmetric representation corresponding 
to the maximum total spin S = N/2. In this subspace 



2 



spanned by the Dicke states [3 |M) = \N/2, M) with 
M = -N/2, • ■ ■ , +N/2, one has: 



ff|M) 



N 



(1 + 7) (iV 2 /4 - M 2 ) - 2M 



\M) 



{ a M-l a M\ M - 2 ) + a i 

A(l- 7 ) 



|M + 2» 



27V 



(5) 



where 



yJ{N/2){N/2 + 1) - M(M ± 1). In the 
following, we will denote by £± the orthogonal sub- 
spaces spanned by the Dicke states |M), such that 
ni^ll-M) = ^l-^O which corresponds to even or odd 
values of (N/2 - M), 

When A < A c and for any 7, the ground state is non- 
degenerate. By contrast, for A > A c , the ground state is 
doubly degenerate in the thermodynamical limit for any 
7 ^ 1 but remains unique in the isotropic case (7 = 1). 
In this limit, the magnetization (per spin) in the z direc- 
tion of the ground state is simply given by : 



N 
1 

N 



1 

2 
1 

2A 



for A < Ar 



for A > A c 



(G) 
(7) 



for all 7. 

To analyze the entanglement properties of the ground 
state we have computed for several values of 7, 
the concurrence introduced by Wootters 0| which, is 
defined as follows. Let us denote by p the reduced 
density matrix obtained from \ip) by tracing out over 
(N — 2) spins. Of course, in our system, the choice of 
the two spins kept is irrelevant contrarily to the ID Ising 
model Next, we introduce the spin- flipped matrix 

p = <Jy®<JyP*<Jy® a y where p* is the complex conjugate 
of p. The concurrence C is then defined by: 



C = max {0, pi — p2 — M3 — A*4} ■ 



(8) 



where the pj are the square roots of the four real eigenval- 
ues, classified in decreasing order, of the non Hermitian 
product matrix pp. This concurrence vanishes for an un- 
entangled two-body state whereas C — 1 for a maximally 
entangled one. Finally, since H couples every spin with 
each other, the two-body entanglement is somewhat "di- 
luted" between all spins, and eventually goes to zero in 
the thermodynamical limit. To get nontrivial informa- 
tions about the entanglement, it is thus crucial to con- 
sider the rescaled concurrence Cr = (N — 1)C where the 
prefactor is simply the coordination number of each spin. 
In symmetric multi-qubit systems, this rescaled concur- 
rence has recently been related to the spin squeezing pa- 
rameter lial 



* 4(Ag g J 2 
N 



(9) 



which measures the spin fluctuations in a correlated 
quantum state ^J. The subscript n± refers to an axis 
perpendicular to the mean spin (S) where the minimal 
value of the variance is obtained. More precisely, for any 
state belonging to £ + or one has: 



e = i - c h 



(10) 



if the matrix elements of the reduced density matrix p 
written in the standard basis {| ft), I ti)>| lt)>| 11)} 
satifies: \pi4\ > P22 [l^l llfljj. In the opposite case, the 
states are not spin squeezed (£ 2 = 1). 

Let us first recall the results in the isotropic case 7 = 1 
which is exactly solvable. As it can be straightforwardly 
obtained from Eq. (|SJ), the (nondegenerate) ground state 
is the Dicke state \N/2) for A < A c , and switches from 
one state \M) to a state \M' < M) as A increases [l7j . 
The concurrence of a Dicke state |M), can be determined 
analytically |2oLl2l| : 



Cr 



— {N 2 - AM 2 
2N 1 



(11) 



y/(N 2 - 4M 2 )[(N - 2) 2 - 4M 2 ]} 



In the thermodynamical limit, this rescaled concurrence 
vanishes for A < A c , jumps to 2 at the critical point A c , 
and decreases, by discrete steps, to 1 at large A . This 
singular behaviour at the transition point is similar to the 
one obtained in the ID Ising model [liEElE]]' except 
that in this latter case, the nearest-neighbour concur- 
rence C(l) is not maximum at the critical point. Con- 
cerning the spin squeezing, its behavior is trivial since all 
Dicke states are not spin squeezed , and no singularity 
can thus be observed on this quantity at the transition. 

For 7^1, the situation is more complex. Indeed, as 
mentionned above, the ground state is doubly degener- 
ate for A > A c in the thermodynamical limit so that, in 
this region, we should, in principle, study the entangle- 
ment obtained from the thermal density matrix (at zero 
temperature) 



Pth. 



(I+X- 



-X-l), 



(12) 



where |+) and |— ) are any two orthogonal ground states. 
Of course, for finite N, the ground state is nondegenerate 
and lies, depending on A, either in £ + or in In the 
thermodynamical limit, the reduced density matrices p± 
built from the corresponding ground state \ip±) by trac- 
ing out over (N — 2) spins become identical. Therefore, 
we have analyzed the entanglement of the true finite N 
ground state (|V'+) or 

We have displayed in Fig. ^ the rescaled concur- 
rence of the ground state, as a function of A for vari- 
ous anisotropy parameters < 7 < 1. For all 7, the 
rescaled concurrence Cr develops a singularity at the 
critical point A c as already pointed out for 7 = 1. How- 
ever, as it can be seen in Fig. |3 the rescaled concurrence 
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FIG. 1: Rescaled concurrence of the ground state as a function FIG. 3: Finite N behavior of 8\Cr for 7 = near the critical 
of A for different values of 7 and for N = 1000. point A c = 1. 



goes to 1 in the thermodynamical limit contrarily to the 
isotropic case where it jumps to 2. More precisely for all 
7^1, one has: 

I-Ck(Am) ~ Ar°- 33±0 - 01 , (13) 
A M -A C ~ 7V-o-66±o.oi 7 (14) 

where Am is the value of A for which C R is maximum. 
In the thermodynamical limit, Cr(Am) goes to 1 while 
Am goes to A c so that the ground state is maximally spin 
squeezed at the critical point (£ 2 = 0). 
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FIG. 2: Behaviors of 1 — Cr(Xm) (upper curves) and Am — 1 
(lower curves) as a function of N for 7 = (o) and 7 = 1/2 
(A). 

To analyze the formation of the singularity at A = A c , 
we have focussed on the case 7 = and plotted in Fig. |3 
the behavior of d\C R near the critical point, for different 
values of N. 

Denoting by A M (respectively A m ) the value of A for 
which d\C R is maximum (respectively minimum), one 
has: 

dxC R {X M ) ~ AT - 33±0 -°\ (15) 



d x C R (\'j ~ -7V°-33±o.o3 5 (16) 
A C -A M ~ iV- 1±0 - 01 , (17) 
X m - A c ~ iV- - 6 ™ 1 , (18) 

for all 7 ^ 1 and at large N. In the thermodynamical 
limit, a real cusp-like singularity is thus observed at the 
quantum critical point. We underline that although we 
are not able to exactly compute the exponent giving the 
large TV behaviors of C R , d\C R (X' M ), and d\C R (X' m ), 
we conjecture that it equals 1/3 . Note that it is also 
the one guessed in Refs. ^| 0] for the scaling of the 
magnetization at the critical point. 

The behaviors of C R and d\C R are notably different 
from those observed in the ID case [||. Indeed, in the 
fully connected system considered here, C R and d\C R are 
extremum at A c whereas in the ID Ising model, d\C(\) is 
the only quantity affected by the transition (C(l) is sur- 
prisingly maximum below the critical point). In addition, 
the scaling behavior of the concurrence and of its deriva- 
tive are different in both models. This simply reflects 
the fact that they do not belong to the same universality 
class as it was already known from the calculations of 
critical exponents. 

In the zero coupling limit (A = 0), the rescaled con- 
currence obviously vanishes since the ground state is in 
this case the fully polarized Dicke state \N/2) and, ac- 
cordingly £ 2 = 1. More interestingly, for 7^0, there 
exists another special value Ao(7) for which C R vanishes. 
For A > Arj(7), the rescaled concurrence is nonzero but 
the ground state is not spin squeezed (£ 2 = 1), whereas 
for A < Ao(7), the spin squeezing is given by (JTUJ. This 
behavior of £ 2 is due to a change in the sign of |px4 1 — P22 
which is always negative above Ao(7). For 7 = 0, such a 
situation never occurs and the ground state is always spin 
squeezed. This is a surprising result since it singularizes 
the case 7 = that, from the phase transition viewpoint, 
belongs to the same universality class as the case 7 7^ 
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FIG. 4: Phase diagram for the ground state spin squeezing 
in the plane (7, Ao). Note that for the isotropic case 7=1, 
the ground state which is a Dicke state is never spin squeezed 
(£ 2 = 1 for any A). 



We have displayed in Fig. U|the "critical line" Ao(7). 
Apart from the very specific case 7 = 1 for which the 
ground state is never spin squeezed, this line is given, 
in the thermodynamical limit, by: Ao = l/\/7 E3- 
We emphasize that this formula is also valid for 7 > 1 
though, in this region, the critical point is readily ob- 
tained by a rescaling of the coupling constant and is given 
by A c = I/7. 

At this stage, we do not completely understand why 
the entanglement properties are so strongly affected by 
the quantum critical point. In particular, the extremiza- 
tion of the (rescaled) concurrence does not seem to be a 
generic characteristic since this phenomenon is not ob- 
served in the ID Ising model in a transverse magnetic 
field [HI, at least for C(l) (13 . Note however that in 
both models, the variation of the concurrence is extremal 
at A c . Although in the present case, we have not exactly 
related the scaling exponent of the entanglement to the 
critical exponents, there may certainly exists some deep 
relations between them which deserves further investiga- 
tions. It would also be interesting to analyze the scaling 
of the Von-Neumann entropy which has been, very re- 
cently, related to the central charge of the conformal the- 
ory associated to the ID quantum spin models Hi). 

Several important issues remains opened. In other sys- 
tems displaying a quantum phase transition, the behav- 
ior of the spin squeezing has never been investigated so 
far. It would be worth determining whether it is always 
minimum at the critical point or not. Indeed, if the con- 
currence is not always maximum at the transition, noth- 
ing prevents the spin squeezing to be minimum as it is 
the case in the present study. Another challenging ques- 
tion concerns the quantum dynamics. For nonstationary 
states, one may wonder how the proximity of a quantum 
critical point influences the time evolution of the entan- 
glement. For a simple initial state fully polarized along 



the field direction, we have already some indications that 
in the fully connected system analyzed here, the rescaled 
concurrence vanishes, at larger times, for A > A c (HI- 
Though we cannot assert that it is a generic situation, it 
is likely that the entanglement of all eigenstates is modi- 
fied at the critical point and consequently, the one of any 
quantum states built from them. Such a study would be 
of primer interest in exactly solvable models. 
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The next-nearest-neighbour concurrence (7(2) is indeed 
maximum at the critical point but its first derivative with 
respect to A vanishes. 



